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=  drop  drag  coefficient 
=  injector  diameter 
=  drop  diameter 
=  turbulence  kinetic  energy 
=  drop  burning  rate 
*  length  of  injector  passage 
=  Ohnesorge  number 
=  pressure 
=  radial  distance 
=  drop  radius 
=  Reynolds  number 
=  time 

=  stream  wise  velocity 
=  Weber  number 
=  streamwise  distance 
-  liquid  volume  fraction 

=  rate  of  dissipation  of  turbulence  kinetic  energy 
=  viscosity 
=  density 
=  surface  tension 
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Subscripts 

c  =  centerline  value 

f  =  liquid-phase  property 

g  =  gas-phase  property 

p  =  drop  property 

o  =  injector  exit  condition 

•»>  «  ambient  condition 

Superscripts 

(-)  -  tune-averaged  property 

(')  «  Favre-averaged  property 
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Intraduciifln 


The  drop  and  spray  combustion  properties  of  the  HAN- 
based  monopropellant  LGP  1845  were  studied.  Drop  burning 
rates  were  measured  with  drops  supported  in  a  combustion  gas 
environment  at  pressures  of  0.2-70  MPa.  Some  internal  gasifi¬ 
cation  of  drops  -  causing  swelling,  partial  bursting,  and  mi¬ 
croexplosions  -  was  observed  throughout  this  region  but  these 
disturbances  decreased  with  increasing  pressure.  Effective  drop 
burning  rates  (including  effects  of  both  surface  gasification  and 
bursting)  were  relatively  constant,  ca.  10  mm/s,  and  were  con¬ 
sistent  with  earlier  strand  burning  rate  measurements  of  gelled 
propellant  Pressure-atomized  combusting  sprays  were  studied 
in  combustion  gas  environments  at  pressures  of  3-9  MPa.  The 
liquid-containing  region  was  significantly  larger  than  earlier 
measurements  of  Birk  and  Reeves,  as  well  as  predictions  based 
on  the  locally  homogeneous-flow  approximation  of  multiphase 
flow  theory.  In  conjunction  with  drop  trajectory  calculations, 
based  on  present  measurements  of  drop  burning  rates,  these 
findings  suggest  significant  effects  of  separated  flow  in  com¬ 
busting  HAN-based  monopropellant  sprays.  . 

Mflmfindatim  ( 


Combusting  monopropellant  sprays  have  applications  to 
regenerative  liquid-propellant  guns,  throttable  thrusters,  and  un¬ 
derwater  propulsion  systems.  The  objective  of  the  present  study 
was  to  experimentally  investigate  aspects  of  monopropellant 
spray  combustion,  seeking  to  extend  earlier  theoretical  results 
obtained  in  this  laboratory1*3  Two  spray  processes  were  con¬ 
sidered,  as  follows:  (1)  the  combustion  properties  of  individual 
drops  supported  in  combustion  gas  environments  at  pressures  of 
0.2-7.0  MPa;  and  (2)  the  structure  of  pressure-atomized  com¬ 
busting  sprays  in  combustion  gas  environments  at  pressures  of 
3-9  MPa.  The  new  measurements  were  used  to  assess  the  im¬ 
portance  of  separated-flow  phenomena  within  pressure-atomized 
combusting  monopropellant  sprays,  i.e.,  effects  of  finite  relative 
velocities  and  transport  rates  between  the  phases.  Similar  to  our 
earlier  work,1'3  the  investigation  was  limited  to  a  hydroxyl-am¬ 
monium  nitrate  (HAN)-based  monopropellant  (LGP  1845) 
which  is  of  interest  for  several  high-pressure  monopropellant 
combustion  systems. 

Individual  drop  burning  rates  are  needed  for  fundamental 
consideration  of  the  properties  of  combusting  monopropellant 
sprays.  Earlier  studies  relevant  to  drop  burning  rates  of  HAN- 
based  monopropellants  have  included  measurements  of  strand 
burning  rates4'6  and  the  burning  rates  of  individual  drops  in 
heated  environments.7'11  McBratney4-5  measured  strand  burn¬ 
ing  rates  of  HAN-based  monopropellants  at  pressures  of  7-100 
MPa.  The  propellant  liquid  was  gelled  with  2  weight  percent 
Kelzan  in  order  to  stabilize  turbulent-like  disturbances  of  the  liq¬ 
uid  surface  that  are  normally  encountered  during  stTand 
combustion  tests  at  high  pressures.  The  strand  burning  rates  of 
gelled  LGP  1845  were  high  (ca.  20  mm/s)  and  the  pressure  de¬ 
pendence  was  relatively  weak  (ca.  p01).  A  frothy  region  was 
observed  at  low  pressures,  where  the  thermal  disturbance  of  the 
combustion  wave  extends  an  appreciable  distance  into  the  un¬ 
burned  propellant,  suggesting  significant  reaction  in  the  con¬ 
densed  phase  for  these  conditions.  While  these  results  are  valu¬ 
able,  however,  the  use  of  a  gelling  agent  raises  questions  con¬ 
cerning  its  influence  on  the  process.  Vosen6  measured  strand 
burning  rates  of  two  ungelled  HAN-based  monopropellant!, 
LGP  1 846  and  a  9. 1  molar  solution  of  HAN  and  water,  at  pres¬ 
sures  of  7-30  MPa.  The  burning  rates  of  both  propellants  were 
very  high,  100-250  mm/s,  and  liquid  surfaces  were  clearly  dis¬ 
turbed,  indicative  of  turbulent-like  instability  of  burning  liquid 
strands  normally  seen  at  high  pressures;  therefore,  these  results 
are  difficult  to  interpret  to  find  the  fundamental  combustion 
properties  of  the  propellants. 

Zhu  and  Law7  studied  the  drop  combustion  properties  of 
LGP  1845  and  other  HAN-based  propellants,  in  combustion 
gases  at  1 170  K  and  1  atm.  The  drops  were  observed  to  beat  up 
with  no  radius  change  at  first,  then  gasify  from  the  surface  for  a 
time  (with  surface  regression  rates  of  ca.  0.2  mm/s),  and  finally 
burst  when  the  drop  diameter  had  decreased  by  roughly  15  per¬ 
cent  Beyer8’9  and  Beyer  and  Teague10  studied  the  combustion 
of  LGP  1846  drops  supported  in  nitrogen  at  temperatures  of 
570-920  K  and  pressures  of  0. 1-8.2  MPa.  These  observations 
yielded  results  similar  to  Zhu  and  Law:7  after  a  beat-up  time  and 
a  period  of  relatively  slow  surface  gasification  (0.2  mm/s  at  730 
K  and  1  MPa)  the  drops  often  burst  -  particularly  the  larger 
drops.  Both  sets  of  drop  experiments  suggest  that  bulk  liquid 
reaction  and  microexplosions  may  be  important  for  combustion 
erf  HAN-based  monopropellants  but  drop  environment  tempera- 
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tunes  were  low  in  comparison  to  the  adiabatic  combustion  tem¬ 
perature  of  the  monopropellant,  (ca.  2150  K);  therefore,  the 
drops  may  not  have  ignited  in  a  manner  representative  of  spray 
combustion. 

Earlier  theoretical  work  in  this  laboratory  addressed  liq¬ 
uid  surface  and  spray  properties  of  combusting  HAN-based 
monopropellants.1'3  Analysis  of  liquid  surface  properties,1*2 
indicated  relatively  high  liquid  surface  temperatures  (in  die  range 
800  -  1050  K  for  pressures  greater  than  10  MPa)  and  unusually 
high  pressures  for  the  liquid  surface  to  reach  its  thermodynamic 
critical  point  (250  MPa  with  an  estimated  uncertainty  of  50  per¬ 
cent).  The  high  surface  temperatures  of  the  liquid  surface  pro¬ 
vides  greater  potential  for  significant  effects  of  chemical  reaction 
in  the  bulk  liquid  than  most  monopropellants,  helping  to  explain 
observations  of  microexplosions  reported  in  Refs.  7-10.  Fur¬ 
thermore,  the  high  critical  combustion  pressure  suggests  that 
spray  combustion  of  HAN-based  monopropellants  involves 
subcritical  combustion  with  a  drop-containing  combusting  spray 
for  most  applications. 

The  earlier  analysis  of  combusting  HAN-based  mono¬ 
propellant  sprays,1*3  was  based  on  the  locally-homogeneous- 
flow  (LHF)  approximation  of  multiphase  flow  theory,  i.e.,  the 
assump,:on  that  velocity  differences  between  the  phases  are 
negligible  at  each  point  in  the  flow;11'13  and  the  thin  laminar 
flamelet  approximation  of  turbulent  premixed  flame  theory,  pro¬ 
posed  by  Bray.14*15  Turbulent  mixing  was  estimated  using  a 
Favre-averaged  k-e  turbulence  model,  with  empirical  constants 
established  from  measurements  in  noncombusting  variable-den¬ 
sity  round  jets,16*17  however,  the  constants  used  are  very  simi¬ 
lar  to  early  proposals  based  on  constant-density  turbulent 
flows.18  The  performance  of  the  analysis  was  evaluated  using 
the  measurements  of  Birk  and  Reeves19  for  pressure  atomized 
combusting  LGP  1846  sprays  at  pressures  of  6-8  MPa.  There 
was  encouraging  agreement  between  predictions  and  measure¬ 
ments,  however,  predictions  were  very  sensitive  to  the  degree  of 
flow  development  at  the  injector  exit  which  was  not  known  very 
well;  therefore,  this  assessment  was  not  definitive.  Later  mea¬ 
surements  of  noncombusting  pressure-atomized  sprays  by  Ruff 
et  al.13  established  the  strong  sensitivity  of  spray  properties  to 
the  degree  of  flow  development  at  the  jet  exit  and  observed  rea¬ 
sonably  good  performance  of  LHF  analysis  in  the  dense-spray 
region  (liquid  volume  fractions  greater  than  0.2)  near  the  injectin’ 
exit  for  atomization  breakup.  However,  these  measurements 
also  disclosed  significant  deficiencies  of  LHF  analysis  for  other 
breakup  regimes  and  in  the  dilute  portion  of  the  spray  —  the  last 
being  in  general  agreement  with  other  recent  evaluations  of  the 
LHF  approach  for  dilute  sprays.11*12 


The  present  investigation  sought  to  extend  past  work 
concerning  both  drop  and  spray  combustion  of  HAN-based 
monopropellants.  Drop  combustion  was  observed  using  an  ap¬ 
proach  similar  to  Beyer10  for  pressures  of  0.2-7  MPa,  however, 
the  drop  environment  more  closely  matched  the  gas  temperature 
of  a  combusting  monopropellant  spray.  Measurements  of  spray 
properties  were  undertaken  seeking  to  confirm  the  measurements 
of  Birk  and  Reeves,19  while  considering  a  broader  range  of  ex¬ 
perimental  conditions  at  pressures  of  3-9  MPa.  The  new  spray 
measurements,  in  conjunction  with  both  LHF  analysis  and  drop 
trajectory  calculations  based  on  the  present  drop  burning  rate 
measurements,  were  used  to  assess  the  importance  of  separated- 
flow  phenomena  for  these  flows. 

Drop  Combustion 
Experimental  Methods 

Apparatus.  Figure  1  is  a  sketch  of  the  drop  combustion 
test  apparatus.  The  supported -drop  technique  was  used  with  the 
drops  exposed  to  gases  in  the  post-flame  region  of  a  premixed 
burner  which  was  operated  within  a  pressure  vessel.  The  pres¬ 
sure  vessel  had  an  inside  diameter  and  length  of  130  and  430 
mm  and  was  fitted  with  two  25  mm  diameter  quartz  windows  so 
that  the  drops  could  be  observed. 

The  premixed  burner  had  a  diameter  of  10  mm  with  a 
stainless  steel  screen  (0.17  mm  diameter  wires,  2000  wires/m, 
square  pattern)  to  help  stabilize  the  flame.  The  gas  flow  rates  of 
the  premixed  burner  were  metered  and  controlled  with  critical 
flow  orifices  and  pressure  regulators.  Burner  operating  times 
were  short,  just  sufficient  to  stabilize  the  premixed  flame  and 
complete  the  drop  combustion  test.  Burner  gas  flows  were 
initiated  and  terminated  with  solenoid  valves  while  the  burner 
was  ignited  with  an  exploding  wire.  The  pressure  rise  of  the 
chamber  (measured  with  a  pressure  transducer)  was  small  in  the 
period  when  the  burner  was  operating,  ca.  5  percent;  therefore, 
the  chamber  pressure  was  set  by  backfilling  it  with  air.  The 
properties  of  the  post-flame  region  of  the  premixed  burner 
roughly  approximated  the  temperatures  of  adiabatic  combustion 
of  tile  monopropellant,  but  contained  significantly  lower  con¬ 
centrations  of  water  vapor,  see  Table  1  for  the  combustion 
product  properties  of  LGP  1845  and  the  burner  gases  (denoted 
burner  1  and  2). 

The  drop  support  assembly  is  illustrated  in  Fig.  2.  The 
drops  were  mounted  on  quartz  fibers,  50-150  pm  in  diameter, 
with  the  bottom  end  of  the  fiber  flame  polished  to  a  bead  of 
somewhat  larger  diameter  to  help  support  the  drop.  The  drop 
was  surrounded  with  a  retractable  shield  to  protect  it  from  tran¬ 
sients  when  the  premixed  flame  was  ignited.  Once  the  premixed 


Table  1.  Combustion  Product  Properties* 


Mixture 

LGP  1845b 

Burner  l6 

Simulent  Gases 

Burner  2d 

Spray 

Temperature  (K) 

2150 

2295 

2230 

2790 

Composition  (%  by  volume^ 
H2O 

69.2 

18.8 

18.0 

19.8 

CO2 

12.9 

9.1 

8.9 

~ 

N2 

17.4 

71.2 

72.6 

38.3 

At 

— 

40.9 

*  Computed  for  10  MPa  using  the  Gordon  and  McBride20  algorithm,  but  effects  of  dissociation  are  small. 
6  Reactant  composition  (%  by  mass):  HAN,  63.2,  TEAN,  20;  and  H20, 16.8. 
c  Major  species  only.  Minor  species  include  CO,  H2,  NO,  OH  and  O2. 

6  Volume  flow  rate  of  burner  gases  (cold)  of  6.28  x  10" 3  m3/s. 
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flame  was  stabilized  the  shield  was  rapidly  retracted  by  fusing  its 
wire  retainer  so  that  the  unbalanced  pressure  force  on  the  shield 
forced  it  to  one  side  of  the  pressure  vessel  where  it  was  stopped 
by  a  rubber  cushion. 

Instrumentation.  Drop  diameter  was  measured  as  a 
function  of  time  using  backlighted  high-speed  motion  picture 
photographs.  The  arrangement  of  the  illuminating  and  camera 
system  is  illustrated  in  Fig.  1.  The  drops  were  backlighted  by  a 
continuous  arc  source,  using  a  condensing  lens  to  direct  the  light 
to  a  diffusion  screen  located  at  one  of  the  windows.  The  pho¬ 
tographs  were  obtained  using  a  high-speed  motion  picture  cam¬ 
era  operating  at  roughly  1000  pictures  per  second  which  incor¬ 
porated  an  internal  timing  marker. 

Results  and  Discussion 

Drop  combustion  at  low  pressures  yielded  very  irregular 
variations  of  drop  diameter  as  a  function  of  time  due  to  bubble 
formation  and  bursting  within  the  drops.  Some  typical  results  at 
low  pressures  are  illustrated  in  Fig.  3.  Drop  diameters  are  plot¬ 
ted  as  a  function  of  time  for  five  tests  at  0.51  MPa  with  initial 

drop  diameters  in  the  range  580-770  |im  and  a  300  jxm  diameter 
bead  on  the  quartz  fiber  to  help  support  the  drop.  The  origin  of 
these  plots  is  somewhat  arbitrary  since  the  motion  of  the  re¬ 
tractable  shield  disturbed  the  premixed  flame  causing  it  to  flap 
for  a  time;  therefore,  the  time  when  the  drop  was  finally  sub¬ 
merged  in  the  post  flame  gases  was  uncontrolled  and  variable. 
The  results  in  Fig.  3  show  swelling  of  the  drop  due  to  bubbles 
in  the  liquid  in  every  case.  The  bubbles  would  periodically 
burst,  carrying  off  some  of  the  liquid,  and  occasionally  the 
bursting  of  a  bubble  was  sufficiently  severe  to  carry  off  all  of  the 
liquid.  At  these  low  pressure  conditions,  internal  reaction  and 
bursting,  with  some  mechanical  removal  of  liquid  caused  by  the 
bursts,  appears  to  be  the  main  mechanism  for  the  reduction  of 
the  drop  diameter. 

The  degree  of  drop  swelling  due  to  the  presence  of  bub¬ 
bles  in  the  bulk  liquid,  and  the  severity  of  drop  bursting,  de¬ 
creased  as  the  pressure  was  increased.  Some  typical  results  at 
higher  pressures  are  illustrated  in  Fig.  4.  Drop  diameter  is  plot¬ 
ted  as  a  function  of  time  for  five  representative  tests  at  2. 1  MPa 
with  initial  drop  diameters  in  the  range  520-680  (lm  and  a  200 

pm  diameter  bead  on  the  quartz  fiber.  As  before,  the  time  origin 
is  arbitrary  due  to  effects  of  initial  disturbances  on  the  premixed 
flame.  All  these  conditions  exhibited  some  degree  of  internal 
bubble  formation,  however,  effects  of  bubbles  bursting  were 
relatively  mild  and  complete  bursting  of  drops  was  not  observed 
at  pressures  of  2. 1  MPa  and  higher. 

Reduced  effects  of  internal  bubbles  at  high  pressures  ap¬ 
pears  to  be  largely  caused  by  increased  gas  density  so  that  a 
given  degree  of  bulk  liquid  reaction  yields  a  lower  volume  of 
gas:  this  reduces  bubble  sizes  and  growth  rates  which  tends  to 
reduce  the  severity  of  bursting  phenomena.  Counter  to  this  is 
the  fact  that  liquid  surface  temperatures  tend  to  increase  with  in¬ 
creasing  pressure  for  the  present  range  of  conditions  (reaching  a 
maximum  at  roughly  25  MPa):1-2  this  is  exnected  to  increase 
rates  of  bulk  liquid  reaction. 

The  time  period  of  drop  swelling,  or  relatively  constant 
drop  diameters,  was  irregular  due  to  uncertainties  concerning  the 
time  when  the  drop  was  submerged  in  the  combustion  gas  envi¬ 
ronment  However,  the  period  when  the  drop  diameter  de¬ 
creased  was  analyzed  to  obtain  effective  drop  burning  rates. 
Plots  of  drop  diameter  as  a  function  of  time  in  the  period  where 
the  drop  diameter  is  decreasing  are  illustrated  for  the  present  test 
conditions  in  Fig.  5.  The  origins  of  these  plots  are  arbitrary 
since  the  data  has  been  plotted  to  overlap  in  tlx;  region  where  the 
drop  diameter  is  decreasing.  In  addition,  conditions  where  the 
drops  burst  completely  at  low  pressures  have  been  excluded. 
Results  at  low  pressures  show  wide  variations  due  to  significant 
effects  of  (nibble  swelling  and  bursting  but  the  diameter  traces 


become  more  regular  and  repeatable  at  high  pressures.  These 
data  were  fitted  to  determine  effective  burning  rates  for  the 
drops,  Kp  =  -drp/dt:  the  fits  are  also  illustrated  in  Fig.  5. 

Present  effective  burning  rates  are  plotted  as  a  function 
of  pressure  in  Fig.  6.  These  results  are  for  drop  diameters  in  the 

range  300-1200  pm  and  include  effects  of  both  internal  reaction 
forming  bubbles  which  buret,  mechanically  removing  some  liq¬ 
uid,  as  well  as  conventional  gasification  at  the  surface  of  the 
drop.  This  combination  of  effects  causes  the  effective  burning 
rate  to  be  highest  at  the  lowest  pressure,  where  bursting  domi¬ 
nates  the  process,  and  then  to  show  relatively  little  change  with 
pressure  over  most  of  the  region  considered  during  present  tests. 
The  strand  burning  results  of  McBratney4’5  and  Vosen6  are  also 
illustrated  in  Fig.  6.  The  present  results  are  a  crude  extension  of 
McBratney’s4-3  measurements  of  gelled  propellants  at  higher 
pressures.  The  results  of  Vosen6  are  much  higher  than  the  rest 
of  the  measurements  due  to  effects  of  liquid  surface  disturbances 
of  burning  liquid  strands  at  high  pressures,  noted  earlier. 

Snrav  Combustion 
Experimental  Methods 

Apparatus.  The  present  spray  combustion  test  apparatus 
was  similar  to  the  arrangement  used  by  Birk  and  Reeves.19  A 
sketch  of  the  apparatus  appears  in  Fig.  7.  The  experiments  were 
conducted  in  the  same  chamber  as  the  drop  combustion  tests. 
The  combustion  environment  was  produced  by  filling  the  cham¬ 
ber  with  a  combustible  mixture  and  then  igniting  it  with  two 
sparks  to  achieve  the  combustion  gas  properties  summarized  in 
Table  1  (denoted  spray).  The  pressure  of  the  spray  tests  was 
adjusted  by  varying  the  initial  pressure  of  the  combustible  gas 
mixture  since  combustion  of  this  gas  approximated  a  constant 
volume  process.  The  combustible  gas  mixture  had  temperatures 
that  were  somewhat  greater  than  adiabatic  combustion  tempera¬ 
ture  of  the  monopropellant. 

The  spray  was  pressure  atomized  using  injectors  having 
exit  diameters  of  0.31,  0.58,  1.08  and  1.17  mm.  The  inlet  of 
tile  injectors  had  baffles,  to  control  any  swirl  in  the  liquid,  and 
smooth  entries,  to  reduce  effects  of  cavitation.  Injectors  having 
length-to-diametcr  ratios  of  2, 17  and  42  were  considered  since 
earlier  work  indicated  that  the  degree  of  flow  development  at  die 
injector  exit  influenced  spray  mixing  properties.2’3’13  The  in¬ 
jectors  were  directed  vertically  upward. 

A  test  was  run  by  placing  a  propellant  sample  (3-4  ml)  in 
the  fuel  delivery  tube  and  filling  the  injector  passage  up  to  its 
exit  A  cap  was  then  placed  over  the  exit  to  prevent  gas  inflow 
when  the  chamber  was  filled  with  the  combustible  gas  mixture 
and  further  pressurized  as  this  gas  burned.  The  propellant  flow 
was  initiated  by  venting  nitrogen  from  an  accumulator  into  the 
fuel  delivery  tube  by  opening  a  solenoid  valve.  Once  the  pres¬ 
sure  of  the  propellant  was  greater  than  the  chamber  pressure,  the 
cap  popped  off  and  the  resulting  propellant  flow  generated  a 
spray  in  the  hot  gas  mixture.  The  process  ended  when  all  the 
propellant  was  consumed.  The  injector  passage  continued  to  be 
purged  by  the  nitrogen  flow  from  the  accumulator  for  a  time  be¬ 
fore  the  accumulator  flow  was  ended. 


Instrumentation.  The  combusting  sprays  were  observed 
using  motion  picture  shadowgraphs  as  illustrated  in  Fig.  7. 
Backlighting  was  provided  by  a  flash  lamp  source,  ca.  1  px 
flash  duration,  which  was  synchronized  with  the  camera; 
therefore,  the  image  of  the  spray  was  effectively  stopped  on  the 
film.  The  shadowgraphs  were  recorded  with  a  16  mm  high¬ 
speed  camera  operating  at  roughly  1000  pictures  per  second, 
using  Tri-X  negative  film.  The  camera  optics  yielded  a  25  mm 
diameter  field  of  view;  therefore,  it  was  necessary  to  adjust  tkw 
position  of  the  injector  to  observe  the  full  length  of  the  liquid 
containing  region. 
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Since  propellant  combustion  does  not  produce  particu¬ 
lates  and  gas  temperatures  are  relatively  uniform  in  monopro¬ 
pellant  spray  flames,  the  boundaries  of  the  spray  were  reason¬ 
ably  well  defined  -  similar  to  past  measurements  of  Birk  and 
Reeves.19  The  films  were  analyzed  to  yield  mean  and  fluctuat¬ 
ing  time  averaged  liquid  volume  fractions,  assigning  dark  zones 
to  unburned  liquid  reactant  and  light  zones  to  gaseous  combus¬ 
tion  products.  For  each  test,  10-30  frames  were  available  for 
analysis  during  the  steady  flow  portion  of  the  spray  combustion 
process.  Separating  dark  and  light  zones  was  somewhat 
subjective;  and  since  the  measurements  correspond  to  line-of- 
sight  projections,  they  are  biased  downstream  and  radially  out¬ 
ward  from  correct  point  measurements  of  mean  and  fluctuating 
liquid  volume  fractions.  Predictions  were  analyzed  to  estimate 
the  line-of-sight  biases,  as  discussed  in  the  next  section. 

Test  Conditions.  Test  conditions  for  the  spray  combus¬ 
tion  tests  are  summarized  in  Table  2.  Most  of  the  injector  flows 
correspond  to  fullv  developed  flow  at  the  injector  exit,  which 
due  to  the  Reynolds  number  of  the  passage  flow  corresponds  to 
turbulent  pipe  flow.  Injection  velocities  were  in  the  range  49-65 
m/s:  these  conditions  correspond  to  the  atomization  break-up 
regime,  i.e.,  a  drop-containing  shear  layer  begins  to  develop  at 
the  liquid  surface  immediately  at  the  injector  exit.12’13 

Theoretical  Methods 

Present  measurements  were  compared  with  the  mono¬ 
propellant  spray  combustion  analysis  developed  earlier  in  this 
laboratory. 1  >3  Drop  trajectory  calculations  were  also  carried  out 
in  order  to  help  assess  effects  of  separated-flow  phenomena. 
Both  analyses  are  described  in  the  following. 

Sorav  Analysis.  The  main  features  of  the  spray  analysis 
will  be  only  briefly  described  in  the  following,  original  sources 
should  be  consulted  for  details. 

The  analysis  involves  use  of  the  LHF  approximation  of 
multiphase  flow  theory11*13  and  the  thin  laminar  flamelet  ap¬ 
proximation  of  premixed  turbulent  flame  theory.14, 15  Turbulent 
mixing  was  treated  using  a  Favre-averaged  turbulence 
model.14*17  This  approach  provides  a  useful  limit  since  both 
multiphase  and  chemical  reaction  phenomena  are  controlled  by 
turbulent  mixing  which  minimized  the  empiricism  needed  for 
predictions,  e.g.,  initial  drop  size  and  velocity  distributions, 
chemical  kinetic  properties,  etc.,  are  not  needed  to  define  the 
problem.  The  main  limitation  of  the  LHF  approximation  is  that 
its  use  generally  tends  to  overestimate  the  rate  of  development  of 
sprays,  particularly  in  dilute-spray  regions  far  from  the  injector 


exit.11*13  However,  Ruff  et  al.13  find  that  the  LHF  approach, 
using  the  present  turbulent-mixing  model,  provided  reasonably 
good  estimates  of  mixing  properties  in  the  near-injector  dense- 
spray  region  of  nonevaporating  pressure-atomized  sprays  in  the 
atomization  breakup  regime  —  conditions  that  are  representative 
of  present  tests. 

The  thin  laminar  flamelet  approximation  implies  that  het¬ 
erogeneous  monopropellant  flames  cover  all  liquid  surfaces. 
Except  for  very  near  the  liquid  surface,  the  liquid  is  at  the  same 
state  as  in  the  injector  while  beyond  the  outer  edge  of  the  thin 
flame  the  gas  has  uniform  properties  equivalent  to  adiabatic 
flame  conditions  noted  in  Table  1.  Under  the  LHF  approxima¬ 
tion,  relative  velocities  between  the  phases  (slip)  are  neglected. 

Other  major  assumptions  of  the  analysis  are  as  follows: 
(1)  steady  (in  the  mean)  axisymmetric  flow  with  no  swirl;  (2) 
low  Mach  numbers  with  negligible  potential  and  kinetic  energy 
changes,  and  negligible  viscous  dissipation;  (3)  boundary-layer 
approximations  apply;  (4)  negligible  effects  of  radiant  energy 
exchange;  (5)  equal  exchange  coefficients  of  all  species  and  heat; 
and  (6)  high  Reynolds  numbers,  so  that  laminar  transport  is 
negligible  in  comparison  to  turbulent  transport  Justification  of 
these  assumptions  is  presented  in  Refs.  1  and  3. 

Under  these  assumptions,  flow  properties  can  be  found 
by  solving  governing  equations  for  conservation  of  mass,  mo¬ 
mentum  and  reaction  progress  variable,  in  conjunction  with  sec¬ 
ond-order  turbulence  model  equations  for  turbulence  kinetic  en¬ 
ergy  and  its  rate  of  dissipation.14-15  The  formulation,  all  em¬ 
pirical  constants  used  in  the  turbulence  model,  and  the  method  of 
solution,  can  be  found  in  Ref.  3. 

The  predictions  were  also  used  to  estimate  potential  ef¬ 
fects  of  line-of-sight  biasing  on  the  measured  distributions  of 
tiquid  volume  fractions  using  a  stochastic  approach  developed 
for  radiation  calculations  in  this  laboratory.21  Knowing  the 
time-averaged  probability  density  function  of  the  reaction 
progress  variable  along  paths  through  the  flow,  the  reaction 
progress  variable  was  simulated  for  a  series  of  statistically-inde- 
pendent  eddies  along  the  path.  Counting  the  presence  of  any 
liquid  in  the  path  as  a  condition  which  would  block  the  light, 
giving  a  dark  image  on  the  film,  yielded  estimates  of  time  aver¬ 
aged  mean  and  fluctuating  liquid  volume  fractions  for  the  path. 
This  procedure  has  not  been  calibrated  using  known  flows, 
however,  it  does  provide  at  least  a  qualitative  indication  of  po¬ 
tential  effects  of  line-of-sight  bias. 


Table  2.  Summary  of  Combusting  Spray  Test  Conditions 


Diameter  L/d  Flow  Amb.  Pres.  Pressure  Inj.  Velocity1* 

(nun)  Type*  (MPa)  Drop  (MPa)  (m/s)  Re6  Ohd  We(* 

Radial  Measurements:  ~  * - 


1.17 

17 

FDF 

3.11 

2.02 

1.17 

17 

FDF 

7.07 

1.79 

1.17 

17 

FDF 

9.10 

1.81 

1.08 

2 

SF 

6.83 

1.93 

0.58 

42 

FDF 

3.19 

2.07 

0.58 

42 

FDF 

6.15 

2.76 

Axial  Measurements: 

0.58 

42 

FDF 

3.22 

2.33 

0.58 

42 

FDF 

6.16 

1.81 

0.58 

42 

FDF 

8.99 

2.14 

0.31 

42 

FDF 

6.51 

3.05 

52.7 

12600 

0.021 

70700 

49.3 

11800 

0.021 

62400 

49.7 

11900 

0.021 

63300 

51.4 

11400 

0.022 

62200 

53.2 

6300 

0.030 

35800 

61.6 

7300 

0.030 

47900 

56.6 

6700 

0.030 

40400 

49.5 

5900 

0.030 

30900 

54.3 

6400 

0.030 

37100 

64.8 

4000 

0.041 

27800 

bUnity  flow  coefficient 

dQh  - 14  /  >/  p(do ,  O- 0.0669  kg/s2 


TDF  =  fully  developed  flow;  SF  =  slug  flow. 
°Re  »  PfUcd/pf,  pf  =  0.0071  kg/ms. 
•Wef^pfUo2*!/®. 


Drop  Trajectory  Analysis.  Direct  assessment  of  the  ap¬ 
proximations  of  the  LHF  approach  for  the  monopropellant 
sprays  was  undertaken  using  drop  trajectory  calculations,  similar 
to  the  approach  used  by  Shearer  et  al 22  and  Mao  et  al 23  for 
nonpre mixed  spray  flames.  These  calculations  were  limited  to 
drops  moving  along  the  axis  of  the  spray.  The  drops  were  as¬ 
sumed  to  be  always  in  contact  with  the  gas  phase  which  was 
taken  to  have  the  properties  summarized  in  Table  1.  Estimates 
of  the  gas  velocity  variation  along  the  axis  were  obtained  from 
the  LHF  analysis. 

Drop  trajectory  calculations  were  limited  to  deterministic 
calculations,  ignoring  effects  of  turbulence/drop  interactions; 
therefore,  mean  gas  velocities  from  the  LHF  analysis  were  used 
in  the  governing  equations  for  drop  motion.  Drops  were  as¬ 
sumed  to  be  surrounded  by  gas  immediately  at  the  injector  exit, 
ignoring  the  all-liquid  core  present  in  these  sprays.13  Effects  of 
drop  heat-up  were  also  ignored:  the  drop  radius  was  assumed 
to  decrease  throughout  the  entire  trajectory  at  10  mm/s  --  based 
on  the  results  of  Fig.  6  for  the  present  test  range.  This  high 
burning  rate  implies  that  the  decomposition  flame  is  located  near 
the  drop  surface,  well  within  the  boundaries  of  the  flow  field 
around  the  drop;  therefore,  gas-phase  properties  used  to  estimate 
drop  drag  were  taken  to  be  ambient  gas  properties  and  effects  of 
forced  convection  or  drop  burning  rates  were  ignored.  Other 
aspects  of  the  analysis  were  similar  to  Refs.  22  and  23:  the  flow 
field  around  the  drop  was  assumed  to  be  quasi-steady  ;  virtual 
mass,  pressure-gradient,  Basset  history  and  gravitational  forces 
were  ignored;  swelling  of  the  drops  was  ignored;  and  drop  drag 
was  estimated  using  the  standard  drag  correlation  for  solid 
spheres. 


Under  these  assumptions,  the  governing  equations  of 
drop  motion  along  the  axis  are  as  follows: 

dxp/dt  =  up 

(1) 

ddp/dt  -  -2  Kp 

(2) 

dup/dt  =  -3pg  Cd  1  Up  -  u  1  (up  -  u)/4 

(3) 

where 

Cd=24  (l+Rep^/frj/Rep,  Rep  51000;  CD=044,  Rep>1000  (4) 

The  initial  condition  is  up  =  u0)  dp  =  dpo  and  xp  =  0  at  t  =  0. 
Equations  (1)  -  (3)  were  integrated  using  a  Runge-Kutta  algo¬ 
rithm. 

Results  and  Discussion 

Due  to  strong  background  lighting,  it  was  not  possible  to 
determine  that  ignition  had  taken  place  from  flame  luminosity. 
Nevertheless,  ignition  was  readily  identified  from  the  pressure 
trace:  inert  liquid  (like  water  or  unignited  monopropellant  spray 
liquid)  caused  the  hot  combustion  products  of  the  premixed  gas 
flame  to  be  quenched  which  resulted  in  a  rapid  reduction  of  the 
chamber  pressure;  in  contrast,  energy  release  from  the  com¬ 
busting  monopropellant  spray  caused  an  increase  of  chamber 
pressure  in  the  period  when  the  propellant  was  flowing.  It  was 
possible  to  consistently  ignite  the  spray  at  pressures  as  low  as 
2.7  MPa,  however,  the  bulk  of  spray  measurements  were  ob¬ 
tained  at  pressures  of  3-9  MPa. 

Measured  and  predicted  time-averaged  liquid  volume 
fractions  along  the  axis,  5fc,  are  plotted  as  a  function  of  nor¬ 
malized  distance  from  the  injector  exit,  x/d,  in  Fig.  8.  Both  pre¬ 
sent  measurements  and  those  of  Birk  and  Reeves19  are  shown 
on  the  plot  Predictions  include  direct  values  of  arc  a*  well  as 
results  allowing  for  line-of-tight  bias,  as  noted  earlier. 


Test  conditions  used  by  Birk  and  Reeves19  were  similar 
to  present  test  conditions,  except  that  injector  L/d  were  in  the 
range  1. 2-2.4  and  the  injector  inlet  was  not  rounded  (see  Lee  et 
al.3  for  a  sketch  of  the  injectors),  and  injector  pressure  drops 
were  1. 5-2.0  times  higher  than  the  present  study.  The  motion 
picture  shadowgraphs  of  both  investigations  were  obtained  in  a 
similar  manner  and  were  analyzed  in  this  laboratory.  Each  set  of 
experimental  results  also  exhibits  a  significant  degree  of  internal 
consistency  and  repeatability  when  plotted  in  the  manner  of  Rg. 
8.  Finally,  pressure  traces  indicated  that  measurements  were 
obtained  for  combusting  sprays  for  both  studies.  Nevertheless, 
present  measurements  exhibit  a  much  longer  liquid-containing 
region  than  those  of  Birk  and  Reeves,19  e.g.,  otfC  =  0.5  at  x/d 
roughly  150  and  25  for  the  two  sets  of  measurements.  Specific 
reasons  for  these  differences  are  not  obvious  since  so  many 
features  of  the  two  studies  were  the  same,  however,  changes  in 
injection  properties  offer  the  most  plausible  explanation.  In  par¬ 
ticular,  the  sharper  injector  inlet  used  by  Birk  and  Reeves19 
could  have  caused  cavitation  in  the  injector  passage  resulting  in  a 
more  finely  atomized  spray  with  a  rapid  rate  of  radial  spread. 
Similarly,  the  injector  used  by  Birk  and  Reeves19  did  not  have  a 
flow  straightener  and  swirl  induced  in  the  injector  flow  passage 
could  have  resulted  in  unusually  high  radial  spread  rates;  al¬ 
though  the  fuel-injection  system  only  involved  rectilinear  motion 
and  doesn't  appear  to  be  fundamentally  prone  to  induce  swirl. 
Finally,  Birk  and  Reeves19  employed  somewhat  higher  injector 
pressure  drops  which  would  be  expected  to  yield  smaller  drop 
sizes  in  the  spray;  nevertheless,  spray  conditions  for  both  inves¬ 
tigations  were  in  the  atomization  breakup  regime  and  the  pres¬ 
sure  drop  increase  doesn't  appear  to  be  sufficient  to  explain  the 
differences  seen  in  Fig.  8  based  on  the  relatively  small  effect  of 
pressure  drop  variations  observed  during  this  investigation.  In 
any  event,  extensive  rechecking  of  measurements  using  the  pre¬ 
sent  injectors  could  not  reproduce  the  results  of  Birk  and 
Reeves.19 

Present  measurements  in  Fig.  8  are  roughly  similar  (in 
terms  of  x/d)  at  all  test  conditions,  with  the  downstream  limit  of 
the  liquid-containing  region  at  x/d  ca.  350.  Since  these  results 
involve  a  range  of  pressures  and  injector  diameters,  tins  behav¬ 
ior  suggests  a  mixing-controlled  process  supporting  the  use  of 
LHF  analysis  —  a  conclusion  reached  in  Ref.  3,  based  on  the 
measurements  of  Birk  and  Reeves.19  Closer  examination  of  tile 
data,  however,  reveals  trends  that  suggest  significant  separated- 
flow  effects.  First  of  all,  results  for  the  0.31  mm  diameter  in¬ 
jector  consistently  exhibit  higher  values  of  afc  at  a  particular  x/d 
than  the  0.58  mm  diameter  injector.  This  is  a  separated- flow 
property  since  drop  diameters  are  not  strongly  affected  by  injec¬ 
tor  diameters  while  drops  of  a  particular  size  must  penetrate  a 
certain  distance  in  order  to  disappear  this  causes  in  a  tendency 
for  penetration  distances,  x,  to  be  constant  for  separated  flows 
rather  than  x/d.11  Another  effect  is  that  (Xfc  at  a  particular  x/d  is 
lower  for  a  chamber  pressure  of  8.99  MPa  than  the  other  pres¬ 
sures  considered  for  the  0.58  mm  diameter  injector  this  behav¬ 
ior  parallels  the  effective  burning  rate  results  of  Fig.  6  where 
drop  burning  rates  at  9  MPa  are  higher  than  for  pressures  in  the 
range  3-6  MPa,  which  are  roughly  the  same.  A  final  effect  is 
that  use  of  long  and  short  L/d  injectors  yielded  roughly  the  «««n^ 
results  while  mixing-controlled  flows  would  result  in  much 
faster  mixing  rates  for  the  long  L/d  injector.13 

Predictions  illustrated  in  Fig.  8  are  for  fully-developed 
flow  at  the  injector  exit,  which  corresponds  to  the  bulk  of  pre¬ 
sent  test  conditions.  Effects  of  ambient  pressure,  injector  diam¬ 
eter  and  injector  Reynolds  number  had  little  effect  on  the  predic¬ 
tions;  therefore,  only  single  lines  are  shown  for  results  with  and 
without  the  line-of-sight  bias  correction.  Comparing  predictions 
with  and  without  the  line-of-sight  bias  correction  indicates  sig¬ 
nificant  effects  of  bias  for  intermediate  values  of  a*;  however, 
predictions  of  the  downstream  end  of  the  liquid-containing  re¬ 
gion  are  not  strongly  influenced  by  the  bias. 
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In  view  of  the  bias  uncertainties,  the  predictions  illus¬ 
trated  in  Fig.  8  are  in  fair  agreement  with  the  measurements  of 
Birk  and  Reeves. 19  This  observation  prompted  earlier  encour¬ 
agement  concerning  the  value  of  the  LHF  and  thin  laminar 
flamelet  approximations  for  analyzing  flows  of  this  type.  How¬ 
ever,  comparison  of  predictions  with  present  measurements  im¬ 
plies  that  use  of  the  LHF  approximation  causes  the  rate  of 
development  of  the  spray  to  be  substantially  overestimated,  in 
agreement  with  most  other  evaluations  of  the  LHF  approxima¬ 
tion  for  sprays.11-13 

Radial  profiles  of  time-averaged  liquid  volume  fractions 
at  various  distances  from  the  injector  are  illustrated  in  Fig.  9.  All 
measurements  shown  in  the  figure  were  obtained  during  the 
present  investigation.  Predictions  shown  on  the  figure  account 
for  line-of-sight  bias  and  are  for  fully-developed  flow  at  the  in¬ 
jector  exit  Similar  to  results  along  the  axis,  predictions  were 
relatively  independent  of  test  conditions  and  only  a  single  line  is 
shown  for  each  streamwise  position.  Results  ignoring  line-of- 
sight  bias  are  narrower  than  the  plots  illustrated  in  Fig.  9,  how¬ 
ever,  the  outer  extent  of  the  liquid-containing  region  is  about  the 
same. 

Similar  to  results  along  the  axis,  the  measured  radial 
profiles  are  crudely  similar  for  all  the  test  conditions  when  plot¬ 
ted  in  the  manner  of  Fig.  9.  In  terms  of  r/x,  the  radial  similarity 
variable  of  turbulent  jets,  the  liquid-containing  region  extends  to 
0.05-0.07,  rather  than  0.15  which  is  the  typical  width  based  on 
scalar  properties  in  turbulent  jets.  Predictions  provide  a  fair  es¬ 
timate  of  flow  widths  near  the  injector  exit  but  progressively  fail 
with  increasing  distance  from  the  injector  exit  --  tending  to  over¬ 
estimate  the  rate  of  development  of  the  flow.  This  behavior  is 
similar  to  other  evaluations  of  the  use  of  the  LHF  approximation 
for  both  nongasifying  and  gasifying  sprays. 

Potential  effects  of  separated  flow  are  examined  directly 
by  the  drop  trajectory  computations  illustrated  in  Fig.  10.  Drop 
velocities  and  diameters  along  the  axis  arc  plotted  as  a  function 
of  distance  from  the  injector  for  an  ambient  pressure  of  10  MPa. 
Predictions  of  velocities  along  the  axis  from  the  LHF  analysis 
are  also  illustrated  on  the  plot,  as  a  reference.  Results  are  shown 

for  initial  drop  diameters  of  10,  20,  100  and  200  pm;  drops 
much  larger  than  200  pm  would  be  subject  to  secondary  breakup 
due  to  excessively-high  drop  Weber  numbers.12  Unlike  LHF 
predictions,  drop  trajectory  calculations  depend  on  the  initial  in¬ 
jector  diameter,  as  noted  earlier  the  results  illustrated  in  Fig.  10 
arc  for  an  injector  diameter  of  1.00  mm. 

The  results  illustrated  in  Fig.  10  clearly  show  significant 
effects  of  separated  flow.  The  LHF  predictions  exhibit  a  decay 
of  velocity  beyond  the  potential-core-like  region  which  is 
roughly  inversely  proportional  to  pressure  -  similar  to  single¬ 
phase  jets.  Due  to  the  small  diameter  of  the  injector,  this  results 
in  a  rapid  deceleration  rate.  Only  the  smallest  drops  (initial  di¬ 
ameters  of  10  pm  or  less)  have  sufficiently  fast  response  to  ap¬ 
proach  the  velocities  of  the  continuous  phase  throughout  most  of 
their  trajectory.  With  increasing  drop  size,  the  drops  progres¬ 
sively  overshoot  the  velocity  of  the  continuous  phase  and  only 
approach  it  again  toward  the  end  of  their  life,  when  they  become 
very  small.  Similarly,  the  drops  pass  beyond  the  end  of  the  liq¬ 
uid-containing  region  estimated  by  the  LHF  analysis  (taken  to  be 
otfc  >  10-4  since  liquid  volume  fraction  never  formally  reaches 
zero  in  the  LHF  analysis  due  to  its  statistical  treatment).  Use  of 
the  drop  burning  rate  estimates  of  Fig.  6,  however,  yields  drop 
trajectories  extending  to  x/d  ca.  300  for  initial  drop  diameters  of 
200  pm.  This  is  comparable  to  present  measurements  of  the 
extent  of  the  liquid-containing  region  suggesting  that  the  results 
of  the  drop  trajectory  calculations  are  at  least  reasonable. 

Taken  together,  the  results  of  Figs.  8-10  suggest  signifi¬ 
cant  effects  of  separated  flow  for  combusting  HAN-based 
monopropellant  sprays  over  the  present  range  of  test  conditions. 


In  view  of  the  relatively  modest  variation  of  burning  rate  with 
pressure  seen  in  Fig.  6,  the  insensitivity  of  drop  drag  properties 
to  pressure,11’12  and  the  relatively  high  critical  combustion 
pressure  of  HAN-based  monopropellants  (ca.  250  MPa1’2)  it  is 
likely  that  separated-flow  phenomena  are  important  for  com¬ 
busting  HAN-based  monopropellant  sprays  for  most  of  their 
range  of  application. 


Conclusions 

The  present  study  involved  measurements  of  the  com¬ 
bustion  properties  of  the  HAN-based  monopropellant  LGP 
1845,  both  as  drops  and  sprays  in  combustion  gas  environments 
at  pressures  of  0.2-9  MPa.  The  spray  measurements,  and  drop 
trajectory  calculations  based  on  the  present  drop  burning  rate 
measurements,  were  used  to  evaluate  earlier  analysis  of  com¬ 
busting  monopropellant  sprays  based  on  the  locally  homoge¬ 
neous  flow  and  thin  laminar  flamelet  approximations,  due  to  Lee 
et  al.1’3  Major  conclusions  of  the  study  are  as  follows: 

1 .  Measurements  yielded  effective  drop  burning  rates  of  ca. 
10  mm/s  for  drop  diameters  of  300-1200  |xm  and  pres¬ 
sures  of  0.2-7  MPa.  The  effective  drop  burning  rate  in¬ 
volved  both  reaction  within  the  bulk  liquid  causing  bub¬ 
ble  formation  and  bursting,  dominating  the  process  at 
low  pressures;  and  conventional  gasification  from  the 
drop  surface,  dominating  the  process  at  high  pressures: 
taken  together,  these  effects  cause  burning  rates  to  be 
relatively  independent  of  pressure  over  the  present  test 
range. 

2 .  Present  measurements  of  drop  burning  rates  at  pressures 
of  0.7-7  MPa  are  generally  consistent  with  earlier  strand 
burning  rate  measurements  of  gelled  propellants  due  to 
McBratney4’5  at  pressures  greater  than  10  MPa. 

3.  Present  measurements  exhibited  a  much  larger  liquid 
containing  region  for  combusting  sprays  at  pressures  of 
3-9  MPa  than  the  earlier  measurements  of  Birk  and 
Reeves19  even  though  test  conditions  and  methods  of 

data  analysis  were  similar,  e.g.,  <XfC  =  0.5  at  x/d  roughly 
150  and  25  for  the  two  sets  of  experiments.  Reasons  for 
these  differences  have  not  been  firmly  established  but 
different  injector  passage  conditions,  possibly  leading  to 
effects  of  cavitation,swirl  and  finer  atomization  for  the 
measurements  of  Ref.  19,  have  been  advanced  as  a  pos¬ 
sible  explanation. 

4.  While  earlier  evaluation  of  analysis  using  the  LHF  and 
thin  laminar  flamelet  approximations  appeared  promising 
based  on  the  measurements  of  Birk  and  Reeves;19  cur¬ 
rent  findings  suggest  that  this  approach  substantially 
overestimates  the  rate  of  development  of  the  flow  which 
is  consistent  with  recent  findings  for  other  pressure-at¬ 
omized  spray  processes.11-13  Separated  flow  phenom¬ 
ena  appear  to  be  important  combusting  for  HAN-based 
monopropellant  sprays  over  much  of  their  range  of  ap¬ 
plication. 
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Fig.  1.  Skadi  of  the  drop  cotnbusuoo  apparatus. 
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DIAMETER  It*  I 


i/d 


Fi|.  !0  Drop  trajectory  properties  along  nit.  Note  dm  afc  -  10*4  at 
i/d  -  220  far  the  LHF  cilceUtion. 


